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We discuss order 1/mf and corrections in inclusive semileptonic decay of a 

B meson. We identify relevant hadronic matrix elements of dimension seven and 
eight and estimate them using the ground-state saturation approximation. Within 
this approach the effects on the integrated rate and on kinematic moments are 
estimated. The overall relative shift in V c b turns out about +0.4% as applied to 
the existing fits. Similar estimates are presented for B — > X s + 7 decays. 
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1. Introduction 

The operator product expansion (OPE) for heavy hadron decays has become the standard 
tool for the evaluation of differential decay rates. While the quality of the OPE may vary 
in different corners of the whole phase space, it has been established that for sufficiently 
inclusive observables the OPE yields an expansion in Aqcd/^q (Q is beauty or charm) 
which, in case of beauty weak decays, converges reasonably well, at least judging from the 
low orders calculated up to now. However, it has been argued that the OPE results in an 
asymptotic series with limitations paralleling those for the perturbative series. In particular, 
this implies a deteriorating behavior at sufficiently high orders; therefore, it is well motivated 
to investigate the higher orders of the power expansion. 

The OPE for inclusive decays yields an expansion in Aqcd/^ with the coefficients which 
are themselves series in a s . Consequently we end up with a double expansion in the two pa- 
rameters. Currently the leading power term, the partonic rate is known to order a^, including 
the differential distributions relevant for the calculation of moments [!]. The coefficients of 
even the first nonperturbative corrections are not completely known to order a s ; the one for 
the chromomagnetic correction is only known at tree level, while the coefficient of the kinetic 
operator [2] can be related to the leading power term by considering the OPE for a moving 
B mesonP] Higher order terms in the OPE are only known at tree level where they can be 



1 This relation is sometimes referred to as a reparametrization invariance. 
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directly constructed by the method discussed in [3J. Certain enhanced terms to order I/tuq 
were calculated for the total decay rate to order 0(a s ) [!]• 

The present paper focuses on the higher orders in the 1 / expansion for b—>c inclusive 
semileptonic decay at tree level. As expected, we observe a proliferation of nonperturba- 
tive expectation values starting from 1/mi. Nevertheless, we identify the set of hadronic 
parameters through order 1/rrvt, all those that appear at tree level. 

Since at higher orders the number of nonperturbative parameters becomes too large, a 
straightforward fit to the data to extract them is not possible. To get around this obstacle, 
we suggest a simple way to estimate the required expectation values. This "ground state 
saturation assumption" evaluates higher-dimension matrix elements in terms of //^, /j,q, p^ D 
and p LS . Using the reasonably well known values for the latter from the data we are in the 
position to numerically assess the higher-order parameters in the power expansion. 

In the next section we describe a general method to calculate the higher orders in the 



1/rriQ expansion at tree level. In Sect. 2.1 we apply this to order 1/mi, including the set of 



the nonperturbative expectation values, and in section 2.2 we do the same for the terms at 
order In Sect. [3] we set up the factorization approximation for the heavy quark matrix 

elements and apply it to evaluate the 5-meson expectation values of dimension seven and 
eight operators appearing at 1 / and 1 /mjj . Using these results we perform a pilot numerical 
study of the corrections at order 1/mf and l/m\ to the hadronic and leptonic moments and 
to \V c b\ in Sect. |4j The impact of higher-order corrections turns out sizable in general; in 
particular, the extracted value of the Darwin term may shift significantly. Nevertheless, the 
value of \V c b\ remains stable. Sect. [5] reports a brief summary of the similar analysis for the 
B^X S + 7 decays, including the analytic expressions for the corrections to the rate and the 
photon energy moments. Here the impact of 1/rnf and terms appears to be insignificant. 
Sect. [6] summarizes our results and presents the main conclusions. Certain aspects of the 
derivation related to the factorization approximation are relegated to Appendices. 

2. Calculational Scheme for at Tree Level 

In this section we give a brief summary of our method to calculate power corrections. It 
follows the calculational scheme used in Ref. f3j. The starting point for the calculation is the 
differential decay rate 

dr= G ^ d ^w, u dL^, (i) 

where the leptonic tensor dL^ v contains all information on the leptons including their phase 
space element. The nontrivial nonperturbative dynamics is encoded in the hadronic tensor 
W^y related, by the optical theorem, to the discontinuity of the time-ordered product of two 
weak currents across the cut. One then starts with a correlator of two hadronic currents 

T ^ iq) = 2Mb! d4xe ~ iqx ( B (^THx)ric(x)m^,m\B(p)) ■ (2) 



Here 



IV = 7^(1 -75) (3) 



is the left-handed weak current and q the momentum transfer to the leptons while P denotes 
the B meson momentum, and 

W„ v = 2lmT lxv . (4) 
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Two-index amplitude T^q) or its discontinuity W^ u (q) can be decomposed into five tensor 
structures with the coefficients given by scalar functions of qo = v ■ q and q 2 : 

W^iq) = -wi(q ,q 2 )g IM/ + w 2 (qo,q 2 )v IM v I/ -iw 3 (q ,q 2 )e lll/p \v p q x (5) 
+w^(q , q 2 )q^q u + w 5 (q , q 2 ){q^v v + v^q u ), 

where v = P/Mb is the four- velocity of the decaying B meson. The structure functions 
W/i and u>5 do not affect decay amplitudes with massless leptons, and the fully differential 
semileptonic rate depending on the charge lepton energy E £ is given by 



d 3 r _ G 2 F \v cb \ 



2 



{2q 2 Wl + [4E £ (qo-E e )-q 2 ]w 2 + 2q 2 {2E e -q )w 3 } . (6) 

The tree-level expansion in 1 /mj is most easily set up by looking at the Feynman diagram 
in Fig. [TJ The double line denotes the Green function of the charm quark propagating in 
the background field of soft gluons in the B meson. We 'rephase' the b fields according to 
b(x) — >e~ imb ( vx,> b(x); this makes the actual &-quark momentum operator to become 

p b = m b v + W. (7) 

The phase factor from b(x) in Eq. Q combines with the background-field c-quark propagator 
[5] l/(ilft — m c ) to yield charm propagator 

to be inserted between the nonrelativistic 'rephased' b fields; this corresponds to the b quark 
momentum Eq. Q. 





iribV+k 



Figure 1: Tree level Feynman diagram for the hadronic tensor in inclusive semileptonic B 
decays 

For semileptonic processes at tree level we only need to multiply ^ by the appropriate 
Dirac matrices Q for the left handed current. A calculation of the OPE series to order l/m b 
requires to expand this expression to n th order in the covariant derivative (iD) according to 

Sbgf = ~7 ~ + -^—(-0)-^— + -r*—{-iH>)— —{-0)— — + ■■■ (9) 
f-m c f-m c f-m c f-m c f-m c f-m c 

The covariant derivatives playing the role of the 'residual' momentum, do not commute in 
general; the above expansion takes care of their ordering. Calculating structure functions 
we take discontinuity of Sbgf; successive terms in the expansion yield higher derivatives of 
S(ml—m 2 +q 2 —2m b qo), together with powers of the spinor factor mi,f—^+m c (m c term can be 
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consistently dropped for purely left-handed weak vertices) producing polynomials in q. The 
cumbersome, for high-order terms, general product of many Dirac matrices is simplified using 
a computer routine. 

Finally, the thus obtained operator expansion should be supplemented by the S-meson 
expectation values of the general form 

b a (iD^)....{%D lln )bp (10) 

where the spinor indices of the b fields are shown explicitly. The field is the full four-component 
QCD field yet redefined by a phase factor. Note that performing the OPE this way yields only 
local operators; however, their expectation values still contain a nontrivial mass dependence 
from the small 'lower' bispinor components and from the explicit subleading terms in the 
equations of motions for b{x) that will be discussed in the following, as well as from the fact 
that the B meson states are the eigenstate of the finite-mass Hamiltonian. 

Selecting the basis for independent expectation values to different orders in l/m& may 
be accomplished in different ways by partially reshuffling the power-suppressed pieces into 
the definition of the lower-dimension expectation values. The most familiar example through 
order ^ adr is the chromomagnetic expectation value /jlq which is usually defined as (b(-aB)b), 
yet often the full Lorentz-scalar operator bia^G^b is used instead; these differ explicitly 
by (/°D+/ 9 is , )/2 ?T1 b- One common strategy assumes that the lower components of the spinors 
are explicitly excluded (while passing to the genuine nonrelativistic b fields, i.e. to the upper 
components corrected by power terms, still may or may not be performed). In this study 
we follow an alternative scheme used in Ref. [3] which may not necessary parallel such a 
separation; it generally keeps all four components of the b field and only reduces the number 
of operators applying equations of motion. While in this scheme the generic operators are 
not expressed through the conventional nonrelativistic b fields, it has an advantage of being 
more compact and more native to the tree-level local OPE calculations. 

The enumeration and evaluation of these matrix elements is conveniently performed in a 
recursive fashion, starting with the operators of the highest dimension, i.e. those with the 
maximal number of covariant derivatives. There one can neglect all 1 /m^-suppressed pieces 
and consider them in the static limit; to facilitate passing on to the next step we, however, 
still use the same full QCD fields. In the static limit one has two possible spinor structures 
of opposite parity [6] : 

(B\b a (iD IMl )....(iD lin )b p \B) = ^W--m» + Mpa^,*"-!** ( n ) 

where ^J- and s\ = ^r7A75^r are the generalizations, to the case of arbitrary velocity 
frame, of the unit and Pauli matrices projected onto the 'upper' components. Temporarily 
introducing an arbitrary restframe velocity vector u„, the standard tensor decomposition 
technique expresses the tensor structures A and B through a minimal set of fundamental 
expectation values at a given order. This can be done by contracting the indices in all possible 
ways. In the following, the matrix elements of this set are referred to as basic parameters for 
a certain order in l/m&. We emphasize that such relations are purely algebraic as long as a 



general parameterization like in Eq. ( 11 ) is adopted; they do not imply any assumption about 
heavy quark limit or subleading correction, or about dynamics - they are just the result of 
rotational symmetry of the B expectation values. Consequently, the same way to relate the 
corresponding tensors to the basic parameters can be used for the expectation values of the 
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operators of any dimension. However, the basic parameters encountered at lower dimension 
generally include higher-order pieces which depend on the choice of the basic set. 

For the operators of dimensions lower than the highest, therefore, we have to include in 
consideration originally all possible 16 bilinear bispinor structures including power-suppressed 



ones when generalizing Eq. ( 11 ). Likewise, we use equations of motions for the b field to relate 
certain structures to higher-dimension expectation values; we do this at the second step, at 
the level of invariants thus reducing directly the set of basic parameters. In particular, 
we successively replace the rightmost and leftmost time derivatives (ivD) acting on b (b) 
according to the QCD equation of motion (see Ref. [TJ, Sect. 2); the same is done to reduce 
certain bispinor structures: 

(iDf + \a, v G^ l-f . 

(ivD)b= — b, — b =— b , 7m7, = <V-^ (12) 

where we have also specified the convention used for c^^] Using various identities for the 
Dirac matrices combined with the second of the above equations one can reduce the number of 



spin structures to one spin-singlet and one spin-triplet, cf. Eq. (11); for the latter we choose 
the dual spin matrix —icr^ u . The spin reduction is less automated than eliminating time 
derivatives of the b fields, nevertheless to the order we are interested in, through operators 
with five derivatives, this is straightforward. The reduction algorithm includes effectively 
eliminating the terms which would mix lower and upper components of the spinors; then 
the leading order of a given expectation value is directly given by the number of covariant 
derivatives in the operator. Spin-triplet operators contain four-dimensional —icr^ u (reduced 
to three-dimensional a mn to the leading order) , and spin-singlet ones are free from any Dirac 
structures. The added convenience of such bases is that the corrections from the presence of 
lower components are automatically 1/m? suppressed. However, we bear in mind that the 
classification itself upon spin properties is valid only to the leading order for a given operator. 
Upon completing the reduction we end up with the expression in terms of basic set at each 
order in Aqcd ; the expectation values at a given order generally add 1 / m^-suppressed terms 
proportional to the basic expectation values of higher orders. 

For the meson expectation values Lorentz covariance of the operators proper is undermined 
by presence of the external 6-states selecting the preferred frame; only rotation invariance re- 
mains. Therefore, the time and the spatial derivatives (and indices in general) are physically 
distinct and should be treated differently. Nevertheless, we write the operators in a su- 
perficially Lorentz-invariant form employing the 6-hadron velocity, as a remnant of the way 
the relations to the basic set have been derived. We denote the projector onto the spatial 
components by 

Then the 'usual' basic parameters, those appearing through order Aq CD can be written as 



2M B f4 = -(B\b v iD p iD a b v \B) IF" 

2M B ii% = \{B\b v [iD p ,iD x ] ( - ia a p) b v \B) IT*'!!^ 



2 Note that this convention differs from the one adopted in a number of earlier papers on the subject. 
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2Mb Pd = ^(B\K 



D p , [iD a ,iD x ]\ b v \B) IT>V 
2Mb pis = \{B\b v [iD p , [iD a ,iD x }}{- ia afi ) b v \B) imi^V 



(13) 



where full QCD b fields and states are assumed. Note that although these parameters may 
differ in higher-order 1/mb pieces from the traditionally used four expectation values, the 
difference for (j% and p? G is only \jm\. 

Performing the calculation along the lines described above and contracting the indices 
obtained in the OPE series Eq. ^ with the tensor structures derived for the expectation 
values generalizing Eq. (11) we obtain the formal l/m^ expansion of the -B-decay structure 



functions, in this case for semileptonic decays. Using the basic parameters in Eq. (13) we 
readily reproduce the known results for the structure functions, lepton spectrum and the total 
decay rate through order l/m\ [8, 9j. 

2.1. Expectation values at 

Passing to higher-order power correction requires first to find independent expectation values 
which can form the set of basic parameters to this order. The starting point are all possible 
index contractions of a matrix element of the form 



(B\b v iD p iD a iD x iD s Tb v \B), 



(14) 



where T is either the unit matrix or —i(7p V . However, not all of them can appear independently; 
T-invariance of the correlator ^ ensures that only Hermitian operators are present. This 
also automatically yields a real differential rate, since the coefficient functions are real. 

Taking this into account, we end up with nine dimension-7 operators, four spin-singlet and 
five spin-triplet: 

b v iD p iD a iD x iD s b v \B) \ (u pa U x5 + n px U aS + n p5 U a) 



2Mb mi = 


(B 


2Mb m 2 = 


(B 


2Mb m 3 = 


(B 


2Mb m 4 = 


(B 


2M B ni5 = 


(B 


2Mb m 6 = 


(B 


2Mb m 7 = 


(B 


2Mb m 8 = 


(B 


2Mb m 9 = 


(B 



b v [iD p ,iD a ] [iD x ,iD s ] b v \B) rPVv A 

K [iD p ,iD a ] [iD x ,iD s ] b v \B) rp A rr 5 
b v [iDp, [»£>„ [iD x ,iD s ]j } b v \B) n ffA n^ 

b v [iDp,iD a ] [iD x ,iD 5 ] ( - %a aP ) b v \B) U ap U^v CT v x 
b v [iDp, iD a ] [iD x , iD s ] ( - UTap) b v \B) n^n^IE' 5 
b v {{iDp,W a }, [iD x ,iD s ]}{-ia a p)b v \B) n^n^ 5 

b v {{iDp,iD a }, [iD x ,iD s ]}{- ia aj3 ) b v \B) n^n QA n^ 

b v iDp, \iD a , [iD x , iDs}} ( - i* a p) b v \B) n^n Aa n CT5 . 



(15) 



The commutator of covariant derivatives equals field strength, and all operators above but 
the first one include as a factor an explicit gluon field operator. Therefore, their expectation 
values without at least one gluon would vanish, whether the hadronic state has zero momen- 
tum or not. Boosting the fo-hadron would not generate per se such expectation values. Only 
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mi can emerge from the boost; it can be identified with average [p 2 ) 2 =p A in the nonrelativis- 
tic picture. However, a subtlety must be born in mind when pursuing such an interpretation, 
for fi 2 as defined in Eqs. (13) contains, along with (p 2 ), an extra piece oc (p A )/m 2 which 



has to be included when checking the effect of the boost through the OPE. This 1/m? piece 



comes both from the presence of the lower components of the b fields in fi 2 in Eqs. (13) and 
from the fact that the upper components proper should be corrected at order 1/m 2 by the 



additional operator 1 + 



(oipy 

8m| 



= 1 + 



9B 



8m 2 b 



to represent the true nonrelativistic fields [7]. 



The corresponding correction most simply can be recovered using two identities 
1 1 . 1 



2M B 



(B\b l0 b\B) 



1, 



this yields 



/4 



2M B 



1 



2m\ 



(B\bb\B) 



(p 4 ) + o 



1 



2m 



2^1 



(16) 



2 ' A 



(17) 

The other operators can also be interpreted in a simple manner. The expectation value 
ttt-2 is proportional to (bE 2 b); likewise is m% proportional to (bB 2 b), where E and B are the 
chromoelectric and chromomagnetic fields. The parameter 1714 is related to (p-roiB) or, by 
the gluon field equations of motion, to a combination of rn^ and (g 2 b (pJ) b) where rot is the 
non-Abelian version of the rotator and Jp the color octet flavor-singlet vector current of light 
quarks. 

The spin-dependent operators have similar interpretation: 77J5 oc (b s- E x E b) and m% oc 
(6 s • B x B b) (different components of the chromoelectric and chromomagnetic fields do not 
commute) while 7717 oc ((p) 2 (s-B) — (s-p)(p-B)), mg oc ((p) 2 (s-B)) and mg is a combination of 
(A(s-B)) and (s-BxB). 



2.2. Expectation values at \/m\ 

The number of possible operators with five covariant derivatives and of their expectation 
values further increases. Requiring T-invariance we get 18 new parameters in total, 7 of 
which are spin singlet and 11 are spin triplet: 



2M B n = 


(B 


biD p (iv ■ 


DfiD p b\B) 


2M B r 2 = 


(B 


biDp (iv ■ 


D)iD p iD a iD a b\B) 


2M B r 3 = 


(B 


biDp (iv ■ 


D)iD a iD p iD a b\B) 


2M B ri = 


(B 


biDp (iv ■ 


D)iD a iD a iD p b\B) 


2M B r 5 = 


(B 


biDpiD p 


(iv ■ D)iD a iD a b\B) 


2M B r 6 = 


(B 


b iD p iD a (iv ■ D) iD a iD p b\B) 


2M B r 7 = 


(B 


b iDp iD a 


(iv ■ D)iD p iD a b\B) 


2M B r 8 = 


(B 


b iDp (iv ■ 


DfiD u (-ia pv )b\B) 


2M B r 9 = 


(B 


b iDp (iv ■ 


D) iD v iDp iD p (—ia pu ) b\B 


2M B r w = 


(B 


biDp (iv ■ 


D) iD p iDp iD u (-ia pu ) b\B 


2M B r n = 


(B 


biDp (iv ■ 


D) iDp iD p iD v (-ia pu ) b\B 


2M B r 12 = 


(B 


b iDp (iv ■ 


D) iDp iD u iD p (-ia pu ) b\B 
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2M B r 13 = 


(B 


biDp (iv • 


D) iDp iD v iD p 


-ia pu )b\B 


2M B r u = 


(B 


biDp (iv ■ 


D)iD p iD p iD u 


'-ia pu )b\B 


2M B r 15 = 


(B 


biDp iD v 


(iv ■ D)iD p iD p 


-ia pu )b\B 


2M B r w = 


(B 


biDp iDp 


(iv ■ D) iD u iD p 


'-ia pu )b\B 


2M B r 17 = 


(B 


biDp iDp 


(iv ■ D)iD p iD u 


-ia pu )b\B 


2M B r ls = 


(B 


b iDp iDp 


(iv ■ D) iD p iD u 


[-ia pu )b\B 



Note that at this dimension we have an odd number of derivatives, therefore at least one 
time derivative (ivD) must be present. This derivative can be commuted through the rest 
until it acts on one of the b quark operators, which yields zero to this order. The remaining 
commutators can be expressed in terms of chromoelectric field, so all these matrix elements are 
necessarily proportional to the gluon field. Unlike the case of an even number of derivatives, 
none of such expectation values can be generated by a B boost alone. 



2.3. OPE to orders and l/m{j 

Using the above set of expectation values and the expansion procedure described earlier in 
this section, we derive all five semileptonic decay structure functions of B mesons through 
order They contain S(m 2 — m 2 +q 2 — 2m b qo) an d its derivatives from first through fifth. 

In view of extremely large number of lengthy terms due to polynomial coefficients in qo an d 
q 2 , to different multiple derivatives of the (5-function for each of the expectation values and 
to algebraic expressions resulting from reduction of various Lorentz structures in the generic 
matrix elements to the basic parameters, we do not attempt to present them explicitly. They 
are handled in Mathematica, with the typical size of the file with their definition around 
0.4 MB. 

The actual inclusive observables we are interested in involve kinematic integrations which 
eliminate the ^-function and its derivatives and yield literal 1/m expansion for the effects 
associated with higher-dimension expectation values. We obtain the particular observable 
using Eq. ([6]) and include the cut on lepton energy if imposed. This is likewise done in 
Mathematica; for the moments with a cut we evaluate the expressions numerically assuming 
certain values for the D = 8 expectation values. Where hadronic mass moments are considered, 
we use the expression for M\ 

M x = ( m b ~m 2 c +q 2 - 2m b q ) + 2(M B - m b )(m b - q ) + (M B - m b ) 2 

= ml + 2(M B -m b )e x + (M B -m b ) 2 ; (19) 

the corresponding factor then multiplies the right hand side of Eq. Q. For integer M\ 
moments the factor is a polynomial in M B — m b with the coefficients that are polynomials in 
qo and q 2 . The quantities m 2 and e x have a meaning of hadron invariant mass squared and 
hadron energy, respectively, at the partonic level in the decay of an isolated b quark. 



3. Estimate of the hadronic expectation values 

As we have seen in the previous section, the independent hadronic matrix elements start to 
proliferate at order 1/mf. This would make their complete extraction from data difficult if 
not impossible. However, not all the corrections may be expected equally important, and 
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neglecting numerically insignificant effects reduces the number of new parameters, possibly 
even to a manageable level. Yet this requires estimating the scale of the emerging expectation 
values well beyond naive dimensional guessing. 

In the present paper we are concerned with the heavy quark expectation values at tree level, 
appearing through dimension 8, and in the following we present a method for estimating them 
based on the saturation by intermediate states. 

Our approach starts with representing the matrix elements of the composite operators of 
interest as a sum over the full set of zero-recoil intermediate single-heavy-quark hadronic 
states. While this is still exact, the approximation scheme we shall employ is to retain in the 
sum the contribution of only the multiplet of the allowed lowest-energy intermediate states. 
The quality of this approximation depends on the degree of saturation of the sum by the lowest 
states. The related accuracy in certain cases can be estimated and the whole approximation 
refined. 



3.1. Intermediate state representation for the relevant matrix elements 

Our goal are forward matrix elements of the form 

(B\biD^iD^---iD^Tb(0)\B), (20) 

where T denotes an arbitrary Dirac matrix. The representation is obtained by splitting the 
full chain iD^iD^ ■ ■ ■ iD^ into A = iD fJil iD fJi2 ■ ■ ■ iD^ and C = iD^ k+1 iD^ k+2 ■ ■ ■ iD^ n . 
The intermediate states representation says that 

(B\bACTb(0)\B) = J^J2( B \b Ab (0)\n) " (n\b CT b{0)\B), (21) 

n 

where we have assumed the B mesons to be static and at rest, \B) = \B(p= (Mb,0))), and 
\n) are the single-6 hadronic states with vanishing spatial momentum. 



We present an OPE-based proof of Eq. (21); a more conventional effective-field-theory 



derivation is given in Appendix A. In either way we introduce a ficticious heavy quark Q 
which will be treated as static (i.e., wiq — > oo and normalization point is much lower than 
Wq), and in the OPE approach consider a correlator of the form 

T AC (qo) = f d 4 x e iqoxo {B\iT{bAQ(x)QCTb(0)} \B); (22) 

note that spatial momentum transfer q has been explicitly set to vanish. 

We shall use the static limit for both b and Q, and hence introduce the 'rephased' fields 
Q(x) = e imQq ° x °Q(x) and likewise for b, and omit tilde in them in what follows. The form of 
the resulting exponent suggests to define u = qQ—m^ + mQ as the natural variable for Tac, 
and 2^j^Tac{^) is assumed to have a heavy quark limit. 

With large ttiq we can perform the OPE for Tac(u) at |w| 3> Aqcd still assuming that 
\u\ <^ rriQ and neglecting thereby all powers of 1/rriQ. In this case the propagator of Q 
becomes static, 

iT{Q(x)Q(0)} = ^V(x) 9(x Q ) P exp U J \ dx ) , (23) 
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and yields 

T AC (u) = (B\bA 1 -C l A^Tb\B), (24) 

— co — ttq — iU 

where 7ro = iZ?o is the time component of the covariant derivative. This representation allows 
immediate expansion of Taq{uj) in a series in 1/oj at large \oj\: 

T AC (co) = -J2 (B\bA ( -^-C ^TblB) . (25) 

k=0 

Alternatively the scattering amplitude can be written through its dispersion relation 



TacH = ^~- de — discT j4c (e) , (26) 

2m / n e—LO+iU 



where we have used the fact that in the static theory the scattering amplitude has only one, 
'physical' cut corresponding to positive to. The discontinuity is given by 

i J d 4 x e iex ° (B\bAQ{x) QCTb(0)\B) (27) 



and amounts to 



discT4 C (e) = ]T i fd 4 x e~^ s e^^ (B\bAQ(0)\n Q ) (n Q \QCTb(0)\B) , (28) 

where the sum runs over the complete set of the intermediate states \uq); their overall spatial 
momentum is denoted by p n and energy by E n . 



The spatial integration over d x and integration over time dxo in Eq. (28) yield (2tt) 5 (p n 



and 2ir S(E n — e), respectively. Therefore only the states with vanishing spatial momentum 
are projected out, and we denote them as |n): 

discT AC (e) = ^ 2m 5(e-E n ) (B\bAQ(0)\n) (n\QCTb{0)\B) . (29) 



Inserting the optical theorem relation (|29|) into the dispersion integral (26) we get 

0)\n) {n\QC 
E n —oj+iO 



rp ( v ^ (B\bAQ(0)\n) (n\QCTb(0)\B) 

T AC ( W ) = ^ ^— — ; , (30) 



and the large-w expansion takes the form 

oo ^ 

T AC {u) = -Y, "Shi E E n (B\bAQ(0)\n) (n\QCTb(0)\B) . (31) 



k=0 



Equating the leading terms in 1/oj of Tac^) i n Eq. (25) and in Eq. (31) we arrive at the 
relation 

{B\b AC^Tb{0)\B) = ^ (B\bAQ(0)\n) ■ (n\Q CTb(0)\B) (32) 



which is the intermediate state representation (21). Note that the projector (l+7o)/2 in the 



left hand side can be omitted since the b field satisfies 6 = 6(l+7o)/2 in the static limit. 
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Considering higher values of k in Eqs. (25) and (31) which describe the subleading in 1/ui 
terms in the asymptotics of Tac(u), we readily generalize the saturation relation (|21[): 



(B\b C ^T^b(0)\B) = Y, (E B -E n ) k (B\bAQ(0)\n) ■ (n\Q C ^T^b(0)\B) . 

n 

(33) 

Thus, each insertion of operator (— ttq) inside a composite operator acts as a factor of the 
intermediate state excitation energy. This is expected, for equation of motion of the static 
quark field Q allows to equate 

id QCb(x) = QTT Cb{x) 

for any color-singlet operator QCb(x). At the same time 

id (n\QCb(x)\B) = —(En — Ms) (n\QCb(x)\B) . 

This reasoning is presented in more detail in Appendix A dedicated to a conventional deriva- 
tion of the saturation relations. 

A couple of comments are in order before closing this subsection. Although we phrased 
consideration for the case of expectation values in B mesons at rest, these assumptions are 
not mandatory. The very same saturation by complete set of heavy quark intermediate states 
of a given spatial momentum holds for matrix elements where initial and final states may 
be different, and may have nonvanishing momenta. They neither have to be the ground 
pseudoscalar states, but with arbitrary spin flavor content. 

Likewise, it is worth noting that even the static approximation for b quarks is actually 
superfluous; the saturation by physical intermediate states relies solely upon large mass of 
their Q quarks. The only modification required for finite mj, is taking care of projector 
introduced by the Q-quark propagator. Using the identity 

1 = ^+75^75 

we arrive at the following generalization: 
(B\bACTb{0)\B) = 

[(B\bAQ(0)\n) ■ {n\Q CTb(0)\B) + (B\b A 75 Q(0)\n) ■ {n\Q C To Tb(0)\B)) . (34) 

n 

The similar relation between 7To and the excitation energy is only modified by the mass shift 
of the finite-mass B meson: 

(B\bA4CTb(0)\B) = Y^[(MB-m b )-{M n -m Q )] k x 

n 

((B\bAQ(0)\n) • (n\Q CTb(0)\B) + (B\b A 75 Q(0)\n) • (n\Q C 75 r 6(0) |B)) . (35) 



3.2. Lowest state saturation ansatz 



The intermediate state representation (21 ) still does not assume any approximation aside from 
the static limit for the b quark, yet it may be used to apply a dynamic QCD approximation. 
The one we employ here uses as an input the l?-meson heavy quark expectation values (20) 
of dimension 5 and 6, which are expressed through (? d and p\ s . 
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All operators with four and more derivatives must have an even number of spatial derivatives 
due to rotational invariance. Thus the operators with four derivatives have either four spatial 
derivatives, or two time and two spatial derivatives. Likewise, the D = 8 operators with 
five derivatives may have four spatial and a single derivative, or two spatial and three time 
derivatives. 

We start with the D = 7 operators with four spatial derivatives, and apply (|32|): 



(B\biDjiD k iDiiD m T b\B) =^2(B\biD j iD k b\n) {n\biDiiD m Tb\B). (36) 

n 

The intermediate states \n) in the sum are either the ground-state multiplet B, B* , or excited 
states with the suitable parity of light degrees of freedom. The ground-state factorization 



approximation assumes that the sum in (36) is to a large extent saturated by the ground 



state spin-symmetry doublet. Hence we retain only the contribution of the ground state 
and discard the contribution of higher excitations. In the case of dimension seven operators 
the result is expressed in terms of the expectation values with two derivatives, i.e. fj% and 
fiQ-, matrix elements involving B* are related to them by spin symmetry. We illustrate 
below the compact result of summation over the multiplet of states. The method we use 
is most economic and turns out particularly transparent when generalizing the ground-state 
approximation. A derivation using the more conventional Lorentz-covariant trace formalism 
is given in Appendix [Bj 

Abstracting first from the heavy quark spin we focus on the indices associated with the light 
degrees of freedom [THj; the corresponding ground state is denoted by |Oo) and its (spinor) 
wavefunction by ^o- The ground-state saturation approximation then reads as 

{UolQiDjiDkiDiiDm Q\Q Q ) = iyi^QiDjiDkQ^) (O |Q iDiiD m Q\VL Q ) (37) 

where the summation over the polarizations of intermediate |fin) is assumed; the spin states 
of the initial and final f^o may be arbitrary. Using rj 

2 2 

{VL \QiD,jiD k Q|Q ) = ^%S jk ^ ~ ^%a jk ^ , where aja k = 5 jk +a jk , (38) 
we get (from now on explicit spinors will be omitted) 

/ 2 \2 2 2 

(n \QiD j iD k iD l iD m Q\n ) = ^-S jk S lm - ^^-(6 jk a lm +a jk S lm ) + (39) 

it lo 

gg (5jm5kl—Sjl5km + 5jrn& k l — 5jia km + aj m 5 k i — (Tji5 km ), 

where spin matrices a act on the spinor indices of the hadron wavefunctions (heavy quark 



field Q can be considered spinless). On the contrary, T in Eq. (36) acts on the spin of real 
b quarks of QCD; for this reason the corresponding cr-matrices are denoted by a® where the 
confusion is possible. 

The matrix elements for B and B* mesons in actual QCD, with arbitrary heavy quark 
spin matrices T®, are directly expressed through those for the f^o-states with spinless heavy 
quarks. To obtain them one would take the trace over the spin indices of both the heavy quark 
and of light degrees of freedom, convoluted with the corresponding meson spin wavefunctions 

M = B + B*a, (40) 



3 Note a difference in the definition of three-dimensional a mn compared to four-dimensional a^v, cf. Eq. dl2p 
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where B and B^. are the B- and i?*-meson fields, respectively. Namely, the generic matrix 
element of an operator b OT® b takes the form 

(M\bOT Q b\M'} = -Tr M ] T Q M't , t = <T 2 Z T a 2 (41) 

if the matrix X describes the spinor part of the matrix element of Q^OQ between the 
states, cf. Eqs. (37)-(39). The charge-conjugate spinor matrix X is evidently S itself for S = l 
and —T, for T, = a k . 

The above rule becomes especially transparent where matrix elements over spinless B 
mesons are considered: then M, M' are unit matrices and one simply takes half the trace over 
spin indices; the structures with an odd overall number of spin matrices vanish (i.e., spin of 
the light cloud must multiply spin of heavy quark in B), whereas products of an even number 
are reduced to a numeric factor by assuming that a® = —a (total angular momentum of B 
meson vanishes!), in agreement with Eqs. (41). 

Using Eqs. p9| ), (41) we get the master equations for spin-singlet and spin-triplet B ex- 
pectation values of D = 7: 

1 -(BfbiDjiDkiDtiDr^lB) = < ^8 jk 8 lm + ^^{8 jm 8 U -5 jl 8 km ) (42) 



2Mb 9 J 36 

1 - /J 2 ^ 2 

{B\biDjiD k iDiiD m a ab b\B) = — ^— (dj k 5i a 5 mb -8 jk 5i b 8 ma + Si m 8 ja 5 kb -Si m 5 jb 5 ka ) + 



2Mb ./«..».««■/ 18 

04)! 

36 



[5jm(filb8ka — dlaSkb) ~ 5ji(5 ka 5 rnb — 5 kb 5 rna ) + 
5kl{5ja$mb — 5j b 5ma) ~ <5fcm(<5 'ja^lb~ & jb^la)] • (43) 



The case of the D = 8 operators with four spatial and one time derivative requires minimal 
modification. Time derivative can occupy the second or fourth position, or stay in the center, 
at the third position. In the former case the evaluation proceeds in the same way, one only 
needs to complement Eq. (38) by the similar relation for Darwin and LS operators: 

3 3 

(Q \QiD j iD iD k Q\Qo) = ~ 6 jk ~ ^jk- (44) 
The two related master equations for D = 8 then become 

i 3 2 3 2 

■ -(B\biD J iD iD k iD l iD m b\B) = -^L5 jk 5 lm + (5 jm 5 kl -5 3l 5 km ) (45) 



2M B j 9 J 36 

-{B\biDjiD iD k iDiiD m a ab b\B) = D G S jk (5i a 5 mb -5i b 5 r 



2Mb / lg 

— 7Tr^8lm(fijaSkb — 8jb$ka) 777, — [^jm(Sl b 5 ka — 5i a 5 kb ) — Sji(5 ka S mb — 5 kb 5ma) + 

lo OD 

5kl(dja,5mb — $jb8ma) ~ S km (8j a 5i b — 5j b 5i a )} . (46) 

For the expectation values with time derivative in the middle position the corresponding 
ground-state contribution vanishes; they appear only due to the 'radially' excited states. 
Therefore, in the ground-state factorization we set these to zero: 

{B\biD j iD k iD Q iD l iD m [a] b\B) ^ 0. 
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expression 


m fc ,GeV 4 




expression 


rrik, GeV 4 




expression 


m fc ,GeV 4 


mi 


1 (a4) 2 


0.11 


m 2 




-0.072 


m 3 


-1(m g ) 2 


-0.082 


m 4 


!(/4) 2 +(/4) 2 


0.39 


m 5 




0.048 


m 6 


1 (/4) 2 


0.082 


mj 


-5„2 ,,2 


-0.42 


m 8 


-8/4/4 


-1.26 


m 9 


( M 2 G ) 2 -f/4/4 


-0.40 



Table 1: Expressions and values for the dimension seven matrix elements 



Finally, we need to consider the expectation values of the form {B\biDjiD^iDi [a] b\B) for 
k = 2,3 which evidently belong to the tower of p? q and p^^g- Likewise, their values could 
be considered as the input describing strong dynamics, along with the latter; yet they have 
not been constrained experimentally. The intermediate states saturating such expectation 
values have opposite parity to the ground state (P-wave states) regardless of number of time 
derivatives. The counterpart of the ground-state saturation approximation here is retaining 
the contribution of the lowest P-wave resonance in the sum; then each power of time derivative 
amounts to the extra power of — e, where e = Mp — M#~0.4 GeV. 

In fact, there are two families of the P-wave excitations of B mesons corresponding to 
spin of light degrees of freedom | or ^. The combinations p^ — pQ, p\> + p\si ■•• receive 
contributions only from the ^-family, whereas the ^-family gives rise 

[TO] (the transition amplitude into the lowest | P-state appears to be suppressed). Therefore, 
it makes sense to consider these two structures separately and approximate 

(B\biDj(-iDo) k+1 iDi b\B) 

(BlbiDji-iDof^iDtajtblB) 

Note that assuming Z\n = £3/2 = e implies p^c^ep^ and —p\s — ^p'g^ * ne nrs * relation seems 
to be satisfied by the preliminary values of p^ and p^ D extracted from experiment. 

3.3. Summary for C?(Aq CD ) and <9(Aq CD ) expectation values 

Combining the above relations we can evaluate all the required nonperturbative parameters 
at order \ jm\ and in terms of a few quantities p\, p^, p^ D , p\ s and e taken as an input. 

Tables [l] and [2] list the resulting expressions for the expectation values (at Z\j<i = £3/2 = e), 
and give the corresponding numerical estimates. For the latter we assume e ~ 0.4 GeV, 
/4 = 0.45GeV 2 , /4 = 0.35 GeV 2 , p% = 0.18 GeV 3 , p\ s = -0.12 GeV 3 . 

The values in the tables are not precision predictions, of course, for a number of reasons. 
First, they depend on p\ , Pq, p^ D , p\ s which are themselves only known with limited accuracy; 
the same holds true for the value of e. This aspect, however, is easy to quantify using the 
expression m Tables [T] and [2} 

Secondly, the estimates are formulated in the infinite mass limit for the b quarks, while 
the parameters of the heavy quark expansion should actually include the full mass depen- 
dence. Generally the finite-mass corrections are governed by the parameter /Uh a dr/2mb and 
can be sizable in b hadrons [11] , up to 15%. However, there is a specific suppression of such 
preasymptotic correction in the ground-state pseudoscalar mesons related to the observed 



_ f-k 2 Pp PLS 1 ^k Pd+Pls \ x (A7\ 

- (^3/2 g + £ l/2 g JOjl (47) 

7 k 2 Pd-Pls , T k 2 Pd+ 2 Pls M8l 

- ~ e 3/2 3 ^ e l/2 3 • V 4 ^ 
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expression 


r fc ,GeV 5 




expression 


r fc ,GeV 5 


r l 


f U 


0.029 




~pIp 3 d 


-0.081 


»"3 




-0.020 


r± 


-\pIpd + \p 2 gpIs + i2 pI 


-0.005 


rs 








r& 




0.029 


ri 








rs 


?Pls 


-0.019 


rg 


-pIpls 


0.054 


rio 


P-gPd 


0.063 


rn 


I (p 2 gPd + pIpIs) - haPls 


0.010 


ri2 


-\{VgPd + pIpIs) - l^GPls 


0.004 


ri3 


\ {-PgPd+pIpIs) + \PgpIs 


-0.046 


ru 


I (p-gP 3 d-pIpIs) + IpgpIs + £ 2 pIs 


0.013 


ri5 








rie 








rn 


?Pls 


-0.019 


ria 









Table 2: Expressions and values for the dimension eight matrix elements 



proximity of these states to the so-called 'BPS' regime [12]. Since we deal here exclusively 
with the pseudoscalar ground state, we expect the finite mass corrections to be substantially 
smaller. 

The major issue is thus the validity of the employed approximation for the matrix elements, 
in particular retaining only the relevant lowest states. The degree to which this ansatz is 
applicable depends on the operator in question, and is expected to deteriorate when the 
number of derivatives (i.e. the operator dimension) increases. 

The dominance of the lowest state and suppression of transitions into highly excited states 
typically holds for the bound states with a smooth potential. In field theory for heavy-light 
mesons this question was studied [13] in two-dimensional QCD, the so called 't Hooft model, 
which is exactly solvable in the limit of a large number of colors. In particular, the ground- 
state expectation values for operators with two spatial derivatives were found to be saturated 
by the first 'P-wave' state to an amazing degree of accuracy (we should remind that there is 
no spin in 1 + 1 dimensions, therefore only one, not two P-wave families). In case this also 
applies to real QCD we would expect a good accuracy of the employed factorization ansatz 
for TO2 and 7715 at order 1/mi and a reasonable one for n, r§, rg, m (and some other related 
combinations) at order l/rri|. All such expectation values are simply given by moments of the 
combinations of two small- velocity structure functions 1^3 /2(e) and Wi/2{e) which are positive, 
which strongly constrains the expectation values. 

The situation is a priori less clear with the operators containing four spatial derivatives. 
In the 't Hooft model the effects related to deviations from the ground-state factorization 
were studied and were found to be nearly saturated, again to a very good degree, by the 
first radial excitation [I3jj^] For some of the expectation values at order they can be 

estimated following the reasoning of Refs. |11} 114]. The nonfactorizable contributions taken 
at face value appear to be about 50% of the ground-state one [15] . 

In reality, the effective interaction in full QCD is singular at short distances due to per- 
turbative physics. Hard gluon corrections lead to a slow decrease of the transitions to the 
highly excited states - yet they are dual to perturbation theory. This is taken care of in 



4 The contribution itself turned out quite significant if normalized literally to (/i^) 2 , apparently since ^ was 

— 4 

anomalously small there lacking the factor of 3, the number of space dimensions. If normalized to A it 
was about 3/4. In actual B mesons fi^. is close to A . 
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the Wilsonian renormalization procedure which is assumed in the kinetic scheme. Prom this 
perspective, one can say that the factorization ansatz yields the expectation values at a low 
normalization point [x < e rac j ~ 0.6 GeV, before the channels to radially excited states open 
up; the excitation energy for such lowest resonance states is probably around 700 MeV. Of 
course, the actual //-dependence at such low scale does not coincide with the one derived from 
perturbation theory; at intermediate excitation energies it must be in some respect dual to 
the perturbative one. 

Keeping this in mind it should be appreciated that even for definitely positive correlators, or 
those expectation values where all intermediate states contribute with the same sign, simply 
adding the first excitation to the ground state contribution may already constitute some 
overshooting. Indeed, a resonance state residing at mass e ra d may be dual to the perturbative 
contribution over the domain of masses 

// 

e rad e rad + ^ rad 

2 2 

where e" ad is the mass of the second excitation, and likewise for higher resonances. Then 
a better approximation for the expectation value normalized at fj, = e ra d would be to add 
only a half of the first excitation contribution. In practical terms, as long as account for the 
power mixing in the perturbative corrections to the conventional Wilson coefficients (most 
notably, of the unit operator) has not been extended to order 1/mf and the effective 

nonfactorizable piece may turn out even less. 

An additional feature of actual QCD is existence of the low-mass continuum contribution 
beyond pure resonances, most notably states like B^tt and their SU(3) siblings. Their 
contribution is typically 1 /N c suppressed and usually does not produce a prominent effect in 
quantities which are finite in the chiral limit (and the corrections to factorization for higher- 
dimension expectation values are). They can be expected to contribute up to 25% of the 
ground state, yet this may be partially offset once the actual QCD conventional expectation 
values are used in the factorization ansatz, that likewise incorporate such states |14j . 

Considering all these arguments, we tentatively assign the uncertainty in the factorization 
estimate to be at the scale of 50%. Yet this should be understood to apply to the 'positive' 
operators where the lowest-state contribution does not vanish and the excited state multiplets 
yield the same-sign contribution. The corrections to factorization will be addressed in more 
detail in the forthcoming paper |15j . 

4. Numerical Estimates for the Higher Order Corrections in the 
Rate and Moments 

Armed with the numerical estimates of all the required expectation values we are in the 
position to evaluate the higher-order power corrections to inclusive B decays. The primary 
quantity of interest is the total semileptonic width r s i(6— >c) used for the precision extraction 
of \V cb \. 
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4.1. T{B^X c iu) 

Assuming the fixed values of and m c we find the following power corrections at different 
orders in l/m^: 

' = -0.043 ' = -0.030 

J- tree tree 

7 = 0.0075 = 0.007 1 = 0.006 , (49) 

J- tree *■ tree i tree 

where r tree includes the phase space suppression factor of approximately 0.63. We have shown 
separately the contribution scaling like l/mfm^. and denoted it by <5r ic . As anticipated [16], 
it dominates the high-order effects and may even exceed the 1 / mf correction, yet it is to some 
extent offset by the regular 1/mjj terms. 



The numerical results (49) suggest that the power series for r s i(6— >c) is well behaved, and 
is under good numerical control provided the nonperturbative expectation values are known. 
Higher order terms induce decreasing corrections except where anticipated on theoretical 
grounds. The estimated overall shift due to higher-order terms 

6T ^ +Sr ^ ^ Q (5Q) 
■L tree 

is well within the interval assessed in |16j and, taken at face value, would yield a 0.65% direct 
reduction in \ V c b\- 

This would not be the whole story, however, for the quark masses determining the partonic 
width are not known beforehand with an accuracy required to extract \V c b\ with the percent 
precision. Rather, their relevant combination is extracted from the fit to the data on kinematic 
moments of the B — > X c tv decay distributions, that in turn are affected by power corrections. 

4.2. Moments 

The key in the OPE evaluation of r s i(&— >c) and, therefore, in extraction of |V^,| are the first 
moments of lepton energy (Ee) and of hadron invariant mass squared (M^), which pinpoint 
the precision value of the combination of and m c that drives the total decay probability. 
Moreover, analyzed through order 1/mjj these two moments turned out to depend on nearly 
the same combination of the heavy quark parameters. This allowed for a nontrivial cross check 
of the OPE-based theory prediction [IT]: essentially, (M|) could be predicted in terms of (Ee) 
and vice versa, once the heavy quark parameters were allowed to vary within theoretically 
acceptable range. 

Therefore, besides the practical question about the shift in the fitted heavy quark param- 
eters, another important issue emerges of whether the consistency between measured values 
of (Ei) and of {M^) persists once higher-order corrections are accounted for. 

Numerically we find 

S(E t ) = 0.013 GeV, 5(M 2 X ) = -0.086 GeV 2 (51) 

where the changes shown are a combined contribution of corrections to order 1 /mf and 1 / . 
In our analysis we follow Ref. [18] and evaluate the moments as the literal ratios 

fdE t dq dg 2 K(E f , q , g 2 ) C(E £ , q , g 2 ) gjg^ 
JdE e d qo d q iC(E e , qo , q i)^£^ 
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without prior expanding the ratio itself in l/m b . Here K is the corresponding kinematic 
observable in question (powers of Eg or Mf-) and C is an explicit kinematic cut if imposed. 
The integrals both in denominator and in numerator are taken directly as obtained in the OPE 
through the corresponding terms in l/m&; in particular, powers of (MB—m b ) for hadronic mass 
moments effectively are not treated as power suppressed. Throughout this paper we perform 
numeric evaluation of higher-order-induced changes in observables discarding perturbative 
corrections altogether; this turned out a good approximation in the kinetic scheme. 

To assess practical significance of the numerical changes in the moments we gauge them 
considering the amount of commensurate shifts required in the values of 'conventional' heavy 
quark parameters used in the OPE so far, to make up for the new effects. Specifically, we 
choose to consider m&, p% and p 3 D for this purpose; changes |<5m&| > 10 MeV, l^ 2 1 >0.1 GeV 2 
and \5p\)\ ^0.1 GeV 3 are deemed significant bearing in mind the estimated accuracy of the 
existing OPE predictions [18]. We do not include m c here (which implies keeping it fixed 
throughout the analysis) for the following reason. The quark mass dependence of the moments 
is essentially given by a combination m b — 0.7m c ; the interval allowed by the fits on individual 
values of mj and m c separately is much wider than 10 MeV, therefore in practice the required 
variation in m c is not independent and is derived from the corresponding variation in m&. 
The stated significance of <5mb~10MeV refers, in fact, to the above combination of masses. 

In actual fits the effect of higher-order corrections is compensated by a change in all heavy 
quark parameters simultaneously. For visualization purposes we, however, quote for each 
moment M the separate values 

SM . 2 6M , 3 5M , KO , 

dm b dfj% dp A D 

as if only one of them were responsible for the adjustment and if the higher-order effect 
had been made up for completely. Clearly, if a particular shift in a heavy quark parameter 
comes out abnormally large, it should simply be discarded: this only signals that the moment 
in question is insensitive to this parameter and the moment rather constrains other OPE 
parameters. On the contrary, if a shift is small, this generally means that the parameter is 



well-constrained and, typically, should not be adjusted. In a sense, the values in Eq. (53) 
would assume that the quality of the fit before including the calculated corrections has been 
perfect, and this clearly is oversimplification. Yet this is suitable to gauge the significance of 
the corrections we study. 

The dependence of the considered moments upon heavy quark parameters, entering de- 



nominators in Eq. ( 53 ) is given in a ready-to- use form in Ref. |18j . From that we obtain 



{E t ): 5m b = —33 MeV (0.022) ; 5/4 = -0.39 GeV 2 (-0.005); Sp 3 D = 0.15 GeV 3 (0.014) 

(Mi): Sm b = -17 MeV (0.011); 5pl = -0.12 GeV 2 (-0.0015); 5p 3 D = 0.086 GeV 3 (0.008) 

(54) 

The dependence of T S \(B) on heavy quark parameters has also been carefully studied [19, 16, 
[18], and above we have supplemented each variation by the corresponding relative shift in 
\V cb \, 



™ * _L ^L^HQP, J_^E!L = ;_ 0.026 GeV" 2 HQP = /i 2 (oo) 

| v^l 2r sl 3HQp * r sl dHQp )_ ni8PoV -3 n J D _5 



0.0013 MeV -1 HQP = m b 
0.026 GeV" 2 HQP = /i 2 
0.18 GeV" 3 HQP = p 3 ) 
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the numbers shown in the parenthesis in Eqs. (54) 



Eqs. (54) suggest that one of the possible 'solutions' is an increase in p^ D by about 0.1 GeV 3 ; 
this value may be affected by a variation in p% of the scale of 0.05 GeV 2 or in mt, about 
lOMeV. The relevance of this solution essentially depends on how precisely the presently 
fitted standard heavy quark parameters accommodate both (Eg) and (M^). 

In actuality, the most precise measurements coming from the threshold production of B- 
mesons at 5-factories require a lower cut on lepton energy. Therefore, the proper analysis of 
the moments should include a E^-cut around 1 GeV. 

Figs. [2]a-f show the size of the nonperturbative OPE terms for the first three central lepton 
energy and hadron mass squared moments depending on the lower cut E*jj nt on the charged 
lepton energy, at different orders in l/m&|^] 



{Ed 




0.010 ■((h i -{h t )'j ) 




(MJk) 





GeV\ 

1.0 1.5 0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 

d e f 

Figure 2: Power corrections to the first three (central) moments of charged lepton energy (upper row) 
or hadron invariant mass squared (lower row), at different orders in the 1/rriQ expansion, 
in units of GeV in the corresponding power. Blue is the effect of an d P% (order 1/m^), 
green at order 1/mf , red at l/fif and magenta finally shows the shift upon including D = 8 
expectation values at order 1/mf 



Keeping in mind that higher moments must generally be sensitive to higher-order OPE 
expectation values, we conclude that at moderate cuts Ef nt < 1.5 GeV preserving sufficient 
'hardness' of the inclusive probability, the power expansion is well behaved; the 1 /m\ effects 
are small compared to the I/ttj^ corrections. This is expected since the IC effects are not 
parametrically enhanced in the higher moments |16} HJ GO] ■ 

At the same time, it is clear that the estimated effects from higher powers in the l/m& 
expansion are not negligible, in particular in the second and higher moments. Those are 
sensitive to p^ and p^ D , and high-order terms may produce their sizable shift. 



To visualize the potential effect, we have plotted the analogies of the shifts in Eqs. (54), 
5mb(Eg ut ), 5fj^(Eg nt ), 5p£ ) (Ef xt ) for all six moments as functions of E^ ut and the correspond- 
ing <5|Vd,|/|Vd,|. Here we present these dependences only for (Ee) and for (Mj^), on single plots, 



separately for 5mb, bpQ and 5p D , Figs. K^a-c; the special role of these two moments has been 



As before, the explicit factor Mb — fUh in hadronic mass moments does not count as a power suppression. 
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Figure 3: Changes in nib, in p% and in respectively required alone to literally offset the effect of 
higher-order power terms in (Ei) (blue) and in (M|-) (magenta), at a given E% ut . In units 
of GeV in the corresponding power 
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(S\V cb \/\Vc\) 
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(0.040) 




(0.014) 






Spi, GeV 2 
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-0.13 


-0.08 


0.33 


(5\V cb \/\V cb \) 


(-0.004) 


(-0.0016) 


(-0.0005) 


(-0.0017) 


(-0.0010) 


(0.0043) 


dpi, GeV 3 


0.16 


0.09 


0.02 


0.09 


0.05 


0.10 




(0.014) 


(0.008) 


(0.020) 


(0.008) 


(0.005) 


(0.009) 



Table 3: Higher order power corrections to the moments with E^ = 1 GeV translated into 
the required conventional heavy quark parameter shifts to offset them; also shown 
are relative shifts in \V c b\ these would induce assuming the fixed value of T S \(B). 
Entries where 5rrib would exceed 100 MeV were left blank 



mentioned earlier in this subsection. The corresponding values (53) for all six moments at a 
representative mild cut E% ut = 1 GeV are shown in Table [3] together with the commensurate 



relative shift <5|Vd,|/|V^,|. The latter is easily estimated using Ref. [16], cf. Eq. (55): 



, - -0.0066 at 5m fe = 10MeV 
1 cbl ~ > 0.0013 at 5/4 = 0.1GeV 2 (56) 



\Vcb\ 



0.009 at 5p 3 D = 0.l GeV 3 



As seen from the plots Figs.[2j the cut-dependence of higher-order corrections shows a gen- 
erally expected behavior which qualitatively follows the behavior already found in the 1 /m 2 
and 1/m 3 effects. The new corrections likewise show mild cut dependence at Ef Ut <lGeV 
and, typically, sharply change above Ey ut ~1.5GeV, in line with the overall deterioration of 
the process hardness with raising cut on E%. 

Based on the numeric pattern of the corrections to the moments we anticipate that inclusion 
of higher-order power-suppressed effects will mostly amount to increase in the fitted value of 
by about 0.1 GeV 3 compared to the fit where only D = 5 and D = 6 nonperturbative 
expectation values are retained, with a possible shift in p\ by about ±0.05 GeV 2 . Figs. @a-f 
illustrate this assertion showing the combined effect of the new power corrections for the 
six moments together with the effect of decreasing p^ D by 0.12 GeV 3 ; the similarity of the 
shifts suggests that the lack of higher-power corrections in the theoretical expressions used 
so far could be to some extent faked by a lower value of the Darwin expectation value. (The 
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third lepton moment is highly sensitive to the Darwin expectation value, and uncalculated 
a s -corrections to the latter may be blamed for the mismatch apparent in the plot Figs. |4]c. 
Besides, lepton energy moments are to a large extent saturated by the parton expressions; 
therefore their high precision allowing to discuss nonperturbative effects in their value relies 
on a high degree of cancellation of conventional perturbative corrections. The extent of such 
cancellation at higher loops is not known beforehand, which warrants a cautious attitude 
towards theoretical precision of higher lepton moments at the required level, and places more 
emphasis in this respect on higher moments of the hadronic mass.) 

A good way to experimentally extract information on the Darwin expectation value is the 
third central hadronic mass squared moment. So far 5-factories did not attempt to measure 
it; it has been extracted with an informative accuracy in the DELPHI analysis |21j . however 
for unspecified reasons this moment was not included into the global fit, according to HFAG. 
We find that the higher-order power corrections (predominantly 1/mt) tend to shift this 
moment by about twice the DELPHI error bar. 



0.014 
0.012 
0.010 


" S{E g ) 




0.006 






0.004 
0.003 


«5<(^-(,e f » 3 > 




0.008 
0.006 






0.004 






0.002 






0.004 






0.002 






0.001 






0.002 
0.000 




E^GeV 


0.000 




Ef u \ GeV 


0.000 




E}"\ GeV 




Figure 4: Effect of including higher-order power corrections (blue) to the first three moments of 
charged lepton energy, upper row, and of hadron invariant mass squared, lower row, and 
effect of decreasing the Darwin expectation value by 0.12 GeV 3 (green), in GeV to the 
corresponding power 



Combining the increase in p^ D with the direct effect on r s i, Eq. (50), we expect an overall 
increase in \ V C \,\ by something like 0.4 percentage points: 



S\Vd>\ 
I Kb I 



+(0.003 -=- 0.005) 



This estimate is based on the expectation values in the ground-state factorization approxi- 
mation and would scale with their magnitude; the actual number may be up to a factor of 
1.5 larger. We emphasize that this would remain only an educated expectation, for the result 
strongly depends on the details of the existing fit to the data, on the precision of different data 
points and on their correlations. The final conclusions should be drawn through incorporating 
the new corrections in the actual fit to the data. 
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5. Power corrections in b— ys + 7 



As another application of the technique described in Sect. [2] we have considered higher-order 
power corrections to the decay rate and to the photon energy moments in radiative b— > s + 7 
decays. They have been treated in the approximation of the local weak vertex; no operators 
with charm quarks or chromomagnetic b—>s vertex were considered. The analysis parallels 
that of the semileptonic case, except that the simplicity of the kinematics (corresponding to 
q 2 = in the latter) leads to reasonably compact analytic expressions even in higher orders. 
We have performed the calculations for an arbitrary mass ratio m s /m;,, however quote here 
the results only at m s = 0. 

Therefore, we assume the b— >s + 7 transition to be mediated by the effective vertex 



A 



^(1-75)6^. 



(57) 



In this approximation power corrections to the integrated decay rate become 
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The first moment, the average photon energy in the decay, is given by 
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For the second and third moments we, to simplify the expressions, quote the corrections to 
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the moments with respect to rrib/2 rather than for the usually considered central moments: 
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The terms through D = 6 in Eqs. (59) and (60) coincide with the known ones, cf. Ref. 



Numerical aspects have been analyzed similarly to the semileptonic moments, employing 
the factorization approximation of Sect.[3]for the expectation values. The corrections turn out 
rather small not only in the integrated width and in the average photon energy, but also in the 
second and even in the third moments corresponding, in the heavy quark limit, to the kinetic 
and Darwin expectation values, respectively. Moreover, accounting for the D = 8 expectation 
values yields small effect compared to the D = 7 one, except for the second moment where 
both are quite suppressed. Direct evaluation results in 



^(B^Xstj) = " °' 036l K " °- 0053 i/ m ? + 0-00064 1/m 4 + 0.00015 1/r , 



5(2E 1 ) = 11 MeV 1/m 2 - 14MeV 1/m 3 +3MeV 1/m 4 + 0.7MeV 1/m s 



12<5((£ 7 -(£ 7 )) 2 > = -0.106GeV 2 3 + 0.002 GeV 2 , 4 + 0.0025 GeV 2 

' b I b ' b 



2A5({E 1 -(E 1 )f)= 0.02 GeV?, 4 - 0.0025 GeV* 5 ; (62) 

'b 'b 

the shifts in the second, third and fourth lines here may be interpreted as an apparent change, 
up to the sign, in m&, /x 2 and — pj-,, respectively. 

The small effect on /i 2 and p^ D , far below the level anticipated in Ref. |22j . is somewhat 
surprising. The above numbers probably are smaller than the corrections due to the non- 
valence four-quark operators of the form bssb appearing at order 0(a s ); they were discussed 
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in Ref. |22j . Sect. 4. It is conceivable that the numerical suppression we obtain is partially 
accidental or is an artifact of the factorization approximation for the expectation values. We 
do not dwell further on this here and plan to look into the issue in the subsequent studies. 

6. Conclusions and Outlook 

In this paper we report a detailed study of higher-order power correction in inclusive weak 
decays of heavy flavor hadrons, focused on the semileptonic B meson decays. The calculations 
existed since the mid 1990s included 1/m? and 1/mj^ corrections. We extended the analysis 
to order 1/m 4 and 1/m 5 ; this is still done at tree level, without computing a s -corrections to 
the power suppressed Wilson coefficients. 

The new calculations contain two elements: deriving the OPE terms proper through the 
order in question, and relating the i?-meson expectation values of the new operators to a set 
of hadronic heavy quark parameters of a given dimension. This has been done in the most 
general setting, by calculating the weak decay structure functions of B mesons through order 
1/m 5 . They are expressed in terms of nine new expectation values at order 1/m 4 and of 
eighteen independent expectation values at order 1/m 5 (more operators appear with 0{a s ) 
corrections, cf. Ref. [?]). Structure functions allow one to calculate any inclusive differential 
distribution incorporating arbitrary lepton kinematic constraints. 

The results for the structure functions are rather lengthy, especially so to order 1/m 5 , and 
we do not quote them; they are generated in a computer program and saved as Mathematica 
definition files. They are used in this form to derive the power corrections to various moments 
of the distributions, both without and with a cut on the charged lepton energy. 

Since the fast growing number of new heavy quark parameters in successive orders in 1/m, 
such an analysis would be of little practical value without means to estimate the emerging 
expectation values with a reasonable confidence. We have presented a formal derivation of 
the saturation representation for the heavy quark matrix elements of operators encountered 
in the tree-level OPE which had been used earlier J3], yet often regarded skeptically. Based 
on this exact representation a ground-state factorization has been formulated and we have 
presented the resulting expressions explicitly for all the encountered I?-meson expectation 
values. 

These two elements of the higher-order analysis allowed us to arrive at meaningful numerical 
estimates and in this way to analyze such important questions as the accuracy of the OPE 
related to the practical truncation of the 1 / expansion and the numerical convergence of 
the OPE series. 

We conclude that the power expansion in inclusive B decays derived from the OPE is 
numerically in a good shape as long as the lower cut on lepton energy remains soft. We 
find that, generally, the 1/m 4 corrections are somewhat smaller than those at 1/m^ and 
typically tend to partially offset the latter; the situation is expected to be different for high 
moments of hadronic invariant mass. The 1/mjj terms in usually considered moments are 
already noticeably smaller and this suggests that the truncation error at this order is largely 
negligible in practice. The notable exception is the effect of 'Intrinsic Charm' (IC) on the 
total decay rate which had been argued [pB] to be a potentially significant effect being driven, 
to higher orders, by an expansion in l/m c rather than in 1/m;,. 

We have also presented the analytic expressions for the corresponding higher-order cor- 
rections in the B — > X s + 7 decays. A pilot evaluation suggested that here the effects are 
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numerically insignificant. A more reliable conclusion, however should include estimates for a 
wider range of hadronic parameters. 

Based on the pattern of the numerical corrections found for various moments important 
in the fit to inclusive semileptonic data, we expect an overall moderate upward shift in the 
extracted value of \ V c b\ about a half percentage point. This comes as an interplay of the direct 
downward shift about 0.65% from the explicit corrections to the semileptonic width and a 
larger upward shift due to the change in the heavy quark parameters extracted from the fit 
to the moments. 

The overall scale of the calculated higher-order power corrections shows that, as a rule, 
they are not negligible at the attained level of precision and are in line with the expectations 
laid down in Ref. [18]. Practically speaking, we expect the corrections to be of the 

same scale as the terms from not yet calculated a s -corrections to the Wilson coefficients of 
the chromomagnetic and of the Darwin operators. 

We expect that the main effect of the estimated power corrections in the fit to the semilep- 
tonic data will be an increase in the Darwin expectation value p^ D by about 0.1 GeV 3 , while 
would not change significantly; neither the main combination of quark masses nib — 0.7m c 
shaping the mass dependence will change essentially. This does not apply to nib or m c sepa- 
rately; the residual dependence on the absolute values is subtle and the changes in the central 
value of nib as large as 50 MeV may not have real significance. 

An increase in the Darwin expectation value would be welcomed theoretically; based on 
the small velocity sum rules and assuming their early saturation we expect 

p 3 D (0.6 GeV) « 0.45 GeV • (0.6 GeV) « 0.15 GeV 3 ; (63) 

evolving this estimate to the routinely used normalization scale 1 GeV we would get a value 
somewhere around 0.25 GeV 3 . 

Extracting the subleading heavy quark parameters from the I^ ut -dependence of the mo- 
ments, a procedure effectively introduced in the fits to the data by assuming a strong correla- 
tion of theoretical uncertainties at different E^ nt for a given moment, looks an unsafe option. 
The cut-dependence of the moments is rather mild until the cut is placed relatively high, and 
there all corrections start to inflate degrading the theoretical accuracy of the OPE predictions. 
This is seen in the power expansion, and the similar behavior is expected from the uncal- 
culated perturbative effects. The correction begin to blow up apparently somewhere near 
Ep ut > 1.65 GeV; to remain on the safe side we suggest to rely on theory at Ef nt < 1.5 GeV. 

Note added. When the present write-up was being finalized, P. Gambino informed us about a 
pilot fit to the semileptonic data including the higher-order corrections to the moments quoted 
here. The preliminary result did not follow our expectation that the bulk of the effects reduces 
to the increase in p 3 -,, although the suppression of the shift in the extracted value of \ V c b\ was 
observed. We note, however, that the reported outcome was largely dependent on the assumed 
strong correlations in the theoretical uncertainties at different Ef nt . More reliable conclusions 
can be drawn once a more general fit analysis is performed. 
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7. Appendices 

A. Conventional derivation of the saturation representation 



Here we present a conventional derivation of the representation of Eqs. (32), (33) as a sum 
over intermediate states, based on the usual field-theoretic language in the framework of the 
effective field theory. We start with the case of operator A without time derivatives: 

{B\bACb(0)\B) = Jd 3 x 5 3 (x) (B\bA(0,x)Cb(0)\B); (A.l) 

time and space coordinates are explicitly shown separately where necessary. A, C typically 
contain derivatives; to avoid related ambiguities one may assume that A acts on the left, on 
the b field. 

Using the property of the infinitely heavy quark fields Q in the sector with a single heavy 
quark 

Q{0,x)Q(0) = 5 3 (x)H^ (A.2) 
(which is nothing but to say they are static) we represent the matrix element as 

(B\bACb(0)\B) = Jd 3 x (B\bAQ(0,x)QCb(0)\B). (A.3) 

In this form we can apply the intermediate state saturation to the product of two heavy quark 
operators bAQ(0,x) and QCb(0): 

(B\bACb(0)\B) = I d 3 x ^2{B\bAQ(0,x)\n) (n\QCb(0)\B) . (A.4) 

J n 

The summation over intermediate states includes both summing over internal degrees of free- 
dom for a particular heavy hadron state (this is traditionally assumed by quantum-mechanical 
summation over states where the center-of-mass motion is factored out), and also integration 
over the momentum p of the hadronic system as a whole: 

2=/t^2> \n) = \n m ,p); (A.5) 

n J \ I states 

below we omit the subscript QM when refer to the conventual states with vanishing total 
spatial momentum ('quantum- mechanical' states). 

Assuming i?-meson is at rest, the first transition matrix element in the r.h.s. of Eq. (A.4) 
has the following dependence on x: 

(B\bAb(0,x)\n,p) = j&(B\bAQ(p)\n); (A.6) 

the integration over d 3 x then yields (2ir) 3 5 3 (p) which removes the integration over the center 
of mass momentum p of the intermediate state | n) , and we arrive at 

(B\bACb(0)\B) = ^ r ^2(B\bAQ(0)\n)(n\QCb(0)\B), (A.7) 

n n 

where \n) include only the states at rest; the factor appears if one uses their standard 
relativistic normalization. 



2(3 



To cover the case of time derivatives we apply the above relation to the following product: 
(B\bA4Cb(0)\B) = -^-^(516^(0)^)^104(76(0)^}, (A.8) 



2M n 

n n 



where time derivative acts on the right. Now, for any local operator O the relation holds 

id (n\O(x )\B) = (E n -E )(n\O(x )\B), (A.9) 

which expresses the fact that the time derivative is given by the commutator with Hamiltonian. 
In the case of heavy quark operators one can additionally apply equation of motion for the 
heavy quark field Q: 

ido{n\QCb(x ,d)\B) = {n\Q K Cb(0)\B) + {n\Q n Cb(0)\B) = {n\Q n Cb(0)\B). (A.10) 



Combining this with Eq. (A.9) and repeating, if necessary, this reduction we obtain 

(n\Q4Cb(0)\B) = (M B -M n ) k {n\QCb(0)\B). (A.ll) 



Therefore, Eq. (A.8) generalizes relation ( A.7| ) in the following way: 



(B\bA4Cb(0)\B) = ^Y,( B \bAQmn){M B -M n ) k {n\QCb{<))\B), (A.12) 



in agreement with the OPE result Eq. (33). 



B. Ground-state contribution in conventional Lorentz-covariant trace formalism 

Here we derive the expressions for the sum over ground-state heavy quark symmetry multiplet 
of meson states, encountered in the factorization approximation for the D = 7 and 8 expecta- 
tion values, in the framework of conventional trace formalism. It was developed in particular 
to describe the heavy quark spin multiplets at different velocities. Since only heavy states 
at rest enter our problem, full Lorentz covariance is superfluous; of four spinor components 
used in conventional formalism only two are independent, corresponding to the spin indices 
of f^o-states considered in Sect. |3.2[ 

With the spin of the heavy quark dynamically decoupled in the heavy quark limit, the B 

(*) 

and B* meson wavefunctions at rest M and can be represented as matrices 

M = v / Ms^i75 MW = v ^^(7e A ), (A.13) 

where we have equated the B and B* masses in the heavy quark limit. One of the indices in 
these matrices corresponds to the heavy quark spin and another to that of the light degrees 
of freedom. Spin symmetry yields the well known trace formula for the matrix elements: 

{H'\biD^iD v T b\H) = -Tr [m h TM h A)j| ht ] (A.14) 

where H and H' are either B or B* , Mjj jj' are the corresponding meson wavefunctions 
and Aji| ht encodes dynamics associated with light degrees of freedom. While a complicated 
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unknown hadronic tensor in the general situation, A|if ht takes a simple form for mesons at 
rest: 



a light 



(toO =1 



components other than spatial vanish. In fact, A f amounts to the corresponding matrix 



3 

light 



6 



to 



liv 1 



(A.15) 



elements of the heavy quark states Qq of Sect. 3.2 



In order to evaluate the ground-state contribution for the matrix element (36) we must 
sum over the states of the (B, B*) spin-symmetry doublet. To do this we employ the relation 
valid for arbitrary R 



2M B 



RM 



(*) 



A 2 2 



(A.16) 



expressing the completeness of the B, B* states in the spin space. Evaluating the right hand 
side of (36) in the ground-state saturation with the trace formula (A. 14) and using relation 
flAlo) we find 

^-(BfbiDjiDkblB) {B\biD l iD m Tb\B)+Y J ^ M ^^ B \ liD 3 iD k^Bl) {B^biD^Y b\B) = 



2M B 



2Mb 



-Tr 



MMA^ ht 



Tr 



MYMA 



light 
Im 



2Mb 



-Tr 



MM { * ] Aj ht j Tr [M^rMAlf 



-Tr 



M 



1+ 7o rMA ii g ht l-7o A h g ht 



light » light 



-Tr \MYMA^A jk 



The product of the two hadronic tensors A llght in (A.15) is given by 

04) 2 



(A.17) 



A light a light 
1V lm IV jk 



9 

2 \2 



9jk9lr, 
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{gjkiVlm + iVjkgim) + 



04) 

36 



■ (gjm9kl — gjl9km + 9jmi& kl — Qjlicr km + *<? 'jm9kl ~ io ' j\9kra) ■ (A. 18) 



In the static limit we therefore obtain for spin-singlet and spin-triplet B expectation values 



1 



2Mb 



(B\biDjiD k iDiiD m b\B) 



04 



2\2 



2 \2 



-9jk9lm + 



04) 

36 



(9jm9kl-9jl9km) (A.19) 



2M B 



{B\biD3iD h iDi%D m (-ia ab ) b\B) 



/4/4 

18 



(9jk9la9mb — 9jk9lb9ma + 



2 \2 



9lm9ja9kb — 9lm9jb9ka) + 



04) 

36 



[9jm{9lb9ka— 9la9kb) ~ 9jl(9ka9mb 9kb9ma ) + 
9kl{9ja9mb — 9jb9ma) ~ 9km(9ja9lb — 9jb9la)], (A. 20) 



reproducing Eqs. (42) and (43). 

For the D = 8 operators with four spatial and one time derivative we follow the same route, 
yet one of the two A llght now describes the operator of the form biDjiDQiD^Yb. Denoting the 
corresponding hadronic tensor of light degrees of freedom by 7*^| ht , we express it through the 
Darwin and Spin-Orbit expectation values: 



7? light 



(A.21) 
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Repeating the same steps we end up with 

-^-(BlbiDjiDoiDkiDtiDrn b\B) = -^Lg jh g lm + {g jm g k l- 9jl9km) (A.22) 

1 - p 3 u 2 

——(B\biDjiD G iD k iDiiD m (-ia ab )b\B) = D 10 G gjk (gi a g m b - gibg-ma) - 
ZlVlB lo 

P'lsP^ PgP'ls 

— Y^dimigjagkb-gjbgka) H — ^ — [gjm(gibgka- gi a gkb) - gji(g ka gmb- gkbg ma ) + 

gki(gjag m b-gjbgma) - gkm(gj a gib-gjbgia)], (A. 23) 



identical to Eqs. (45) and (46). 
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